Vacancy-induced magnetism in Sn02: A density functional study 
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We study the magnetic and electronic properties of defects in Sn02 using pseudopotential and all 
electron methods. Our calculations show that bulk Sn02 is non-magnetic, but it shows magnetism 
with a magnetic moment around 4.00 hb due to Sn vacancy (Vsn)- The magnetic moment comes 
mainly from O atoms surrounding Vsn and Sn atoms, which couple antiferromagnetically with the 
O atoms in the presence of Vsn. The coupling between different Sn vacancies is also studied and 
we find that these defects not only couple ferromagnetically but also antiferromagnetically and 
ferrimagnetically. Our calculations demonstrate that the experimentally observed giant magnetic 
moment of transition metal doped Sn02 can be attributed to Vsn. 

PACS numbers: 75.50.Pp, 61.72.Ji, 71.22.+i 



In many diluted magnetic semiconductor (DMS) the 
non-magnetic matrix is a conventional compound semi- 
conductor such as GaAs^ or InAs^. These DMSs have 
low solubility limit and their Curie temperatures (Ic)s 
are well below room temperature (RT), which disqual- 
ify them for spintronic devices. In other classes of DMS 
the transition metal (TM) is embedded in oxide semi- 
conductors, which are conventionally known as oxide- 
DMS (ODMS), such as ZnO with Co or Mn doping^'^-l, 
Ti02 (anatase) with Co*' and Sn02 with Cc^. These 
ODMSs have large magnetic moments and their Tqs are 
well above RT. Therefore these are good candidates for 
spintronic devices. 

Sn02 is a wide band gap material which has been used 
as a transparent conducting electrode in solar cells^ and 
flat-panel display^,. This compound has a rutile struc- 
ture with distorted octahedral coordination. Ogale et alJ 
found that Co- doped Sn02 not only exhibits ferromag- 
netism above RT but also shows a giant magnetic mo- 
ment (GMM) of 7.5±0.5 /is/Co, which was reported to 
be the first GMM observed in any DMS. The GMM was 
attributed to either the cobalt orbital magnetic moment, 
which was probably not quenched, or the appearance of 
some moment on the atoms surrounding the cobalt in the 
matrix. Ferromagnetism with high Tqs was also observed 
in Feia, Cr^, V-i^, and Ni^^ doped-Sn02. The ferromag- 
netism in V, Cr, and Ni doped-Sn02 was found to depend 
on the nature of the substrate on which the samples were 
depositediiii^iii^. More recentlyl^, it has been observed 
that Mn-doped Sn02 also showed very large magnetic 
moment i.e. 20.0 /is at low doping concentration which 
is far above the spin magnetic moment of Mn. 

To date, there are few cases where defects are believed 
to be the origin of ferromagnetism and each case has a dif- 
ferent physical origin. These cases include CaO-^, SiC^^, 
GaBg^, and IIf02^^. Each compound has different crys- 
tallographic environment and local symmetry which gives 
birth to different origins of ferromagnetism in these ma- 
terials. Experimentally it has been observed these ma- 
terials show ferromagnetism without doping of TM, like 
for example IIf02 and CaBg. Our case is quite different 



than those mentioned above. Sn02 has a rutile structure 
in which each Sn is surrounded by an oxygen distorted 
octahedron. Also, Sn02 shows very large magnetic mo- 
ments 20 hb, for Mn-doped Sn02) when doped with 
TM. However, we will show that it is possible to induce 
magnetism in Sn02 without doping of TM. Most impor- 
tantly, the previous studied materials did not consider 
the magnetic coupling between the defects. We study in 
detail that coupling to prove that Sn02 has magnetism 
without TM doping. 

Considering the defect chemistry of Sn02 , it has been 
concluded that the TM cations are in a -1-3 ionic state, i.e. 
Mn+3, Fe+3 and Co+^ii. These TMs have lower valence 
than Sn for which they substitute and therefore will influ- 
ence the defect chemistry of Sn02. On the other hand, 
in structures where ionic or metallic binding predomi- 
nates, as many as half of the cation sites may be vacant 
or may occupy interstitial sitesiS. In other words we can 
say that vacancies are quite common in many crystals, 
both in closed-packed metallic structures as well as in 
open covalent structures^. Although TM doped-Sn02 
shows RT ferromagnetism with large magnetic moments, 
which is ideal for spintronic devices, the origin of GMM 
is still unclear. To flnd the physical origin of GMM we 
considered two kinds of defects i.e. Sn vacancies (Vsn) 
and oxygen vacancies (Vo)- Interestingly we found that 
Vsn has a magnetic ground state with a large magnetic 
moment, which indicates that Vsn can be responsible for 
GMM. On the other hand Vo does not show any mag- 
netism. Detailed energetics and other kinds of defects 
are outside the scope of the present work, but they have 
been studied in detaiP^ before the discovery of ferromag- 
netism in Sn02. So here we will focus only on Vsn and 
Vo in connection with magnetism. 

We performed density functional theory calculations 
using the SIESTA code^^ with a double-zeta polarized 
(DZP) basis set. For the exchange and correlation po- 
tential, we used both the local density approximation 
(LDA)^'^ and the generalized gradient approximation 
(GGA)^. The Hamiltonian matrix elements were calcu- 
lated on a real space grid defined with a plane-wave en- 
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FIG. 1; Total DOS for O-defect (Vo) system calculated with 
LDA. The inset shows the total DOS for a pure Sni6032 super 
ceU. 



ergy cutoff of 300 Ry. A 4x4 x 6 Monkhorst-Pack (MP) 
mesh for one unit cell (6 atoms) and a 3 x 3 x 3 MP for a 
2x2x2 supercell (48) atoms were used for k point sam- 
pling. The atomic positions were relaxed until all forces 
were smaller than 0.05 eV/A. To confirm our results, we 
also carried out calculations using the full potential lin- 
earized augmented plane wave (FLAPW)^ with GGA^i, 
48 special k-points and about 2450 basis functions, which 
yielded qualitatively the same results. As a starting point 
we used our optimized LDA (GGA) parameters of Sn02 
which are: a = 4.715(4.780) A, c = 3.248(3.268) A, and 
u — 0.307(0.307). For Vsn we considered a defect con- 
centration of 0.0625, which corresponded to a single Sn 
vacancy in the unit cell of 48. For Vq we removed one 
O atom from the same unit cell. 

In a perfect Sn02 crystal the Sn is in a 4-1- state due to 
the donation of its 4 valence electrons to the O2 complex, 
which is in a 4- state. This results in Sn02 as a non- 
magnetic and wide band gap material where all bands 
are occupied. The total DOS of pure Sn02 calculated 
with and without spin polarization gives the same result 
and is shown in the inset of Fig. [TJ The valence band is 
formed mainly by O 2p orbitals and is full whereas the 
conduction band is formed mainly by the Sn 5s orbitals 
and is empty. 

On a simple ionic picture, removal of neutral oxygen 
would lead to the reduction of tin from Sn'^^ to Sn^+ 
state, as has been confirmed by earlier spectroscopic 
studies^. Dangling bonds cannot be expected since Sn 
can also exist in a Sn^+ state and it will complete its 
bonding with the nearby O atoms, which will result in 
SnO which is an insulator—. We find that Vq is a non- 
magnetic insulator. The total DOS for Vq is shown in 
Fig. [TJ The oxygen defect creates a defect band inside 
the band gap of pure Sn02 which does not destroy the 
insulating behavior of Sn02, but decreases its band gap. 
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FIG. 2: (a) Total DOS for Vsn in the PM state, (b) Total 
DOS for Vsn in the FM state. The inset shows the DOS 
calculated with FLAPW. Partial DOS of (c) O and (d) Sn 
atoms for Vsn- Dotted and dashed lines show s and p orbitals, 
respectively, whereas the solid lines show the total DOS. 



This indicates that the O defect is electronic in nature. 

When one Sn vacancy is introduced in the supercell 
the bond length of Sn-O is decreased by ~0.14 A after 
relaxation. In this case holes are created at the O sites. 
The DOS for Vsn in the paramagnetic case using GGA 
is shown in Fig. [2] (a). 

We find that the band structure is identical to pure 
Sn02, except close to the Fermi energy. Now, a hole 
band formed by the O surrounding Vs„ is created on 
top of the valence band. Therefore a cation vacancy 
removes symmetrized orbitals from the bonding bands 
and raises them above the Fermi energy. The DOS indi- 
cates that the electrons are partially localized, compared 
to pure Sn02 where electrons are completely localized 
due to strong bonding, and there is a metallic behavior. 
The finite density at the Fermi level indicates that these 
holes will mediate a magnetic interaction, if Vsn is show- 
ing any magnetism. To see more deeply whether they 
correspond to spin up or spin down electrons we spin po- 
larized the band structure. All our calculations (SIESTA 
and FLAPW) show that Vsn is more stable in the mag- 
netic (M) state than in the paramagnetic (PM), with an 
energy difference, £^m — -E-pm, which ranges between -0.53 
(LDA) and -0.71 eV (GGA). 

The GGA spin polarized DOS is shown in Fig. [H (b). 
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FIG. 3: Spin density contours for Vsn in the (110) plane cal- 
culated with GGA using FLAPW. The lowest contour starts 
from 2. Ox 10~* electrons/(a.u)^ and the subsequent lines dif- 
fer by a factor of 2. 
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All the states below -10 eV are spin degenerate and for 
simplicity are not shown. However, the situation changes 
dramatically around the Fermi level. We see the DOS is 
spin split, which implies that the Sn defect induces mag- 
netism in this system. We also find a large magnetic 
moment of 3.85 (4.00) /is with LDA (GGA). The Sn de- 
fect creates holes at the neighboring O ligands which are 
localized for spin up while they are partially localized 
for spin down electrons. So when we introduce one Sn 
defect we are creating two ligand O holes which in turn 
give a total magnetic moment 4.00 /is, as shown by our 
ab inito calculations. We then analyze the PDOS and 
we find that the O atoms surrounding the Sn defects are 
the main ones that contribute to the magnetism. The 
PDOS of one of these O atoms is shown in Fig. [5] (c) . We 
see the spin down state is partially occupied, which con- 
tributes to the magnetic moment of Vsn- For comparison 
we also show the PDOS of a near Sn atom (Fig. [21(d)), 
which produces almost no contribution to the magnetic 
moment. So the magnetic moment is mainly attributed 
to the O atoms due to the dangling bonds that create 
unpaired electrons. The total DOS also shows that ma- 
jority spins keep the insulating behavior of pure Sn02 
whereas minority spins are metallic. Such DOS indicates 
half metallic behavior, which is very important for spin- 
tronic devices. The LDA DOS, which is not shown here, 
produces a similar trend. The FLAPW also shows that 
the ground state of Vsn is magnetic with a magnetic mo- 
ment of 4.00 fiB- For comparison the total DOS calcu- 
lated with FLAPW within GGA is shown in the inset of 
Fig. m (b) . To see the coupling between O and Sn atoms 
for the Vsn system, we calculated the spatial projection 
of the spin densities. The spin contours are shown in 
Fig.H 

We see that the Sn atoms are negatively polarized and 
couple antiferromagnetically with the O atoms and the 
spin density clearly indicates that the magnetic moments 
are localized at O sites. These facts are also confirmed 
by the MuUiken charge analysis (SIESTA) as well as by 
the magnetic moment within the MT spheres (FLAPW). 
The spin projected densities also show that the charge 



FIG. 4: The exchange coupling constant J{eV) as a function 
of Vsn separation. The dots are data points whereas the solid 
line is the RKKY model fitted with fcp = 1.19 A"^ 



of O atoms near Vsn spreads towards the defect due to 
the dangling bonds and the magnetic moments of such 
O atoms are increased compared to those far away from 
Vsn- We can also see that the the magnetic moments 
of O atoms decrease rapidly as the distance of O atoms 
from Vsn increases. 

Until now we focused only on single Sn vacancies 
and we concluded that Vsn is magnetic using FLAPW 
and SIESTA. However, to specify wether there is ferro- 
magnetism or antiferromagnetism, which are cooperative 
phenomena, there should be at least two entities in the 
unit cell that can determine the true magnetic ground 
state properties. We therefore move to discuss the cou- 
pling between the Vsn defects. We considered two Vsn 
defects in the supercell and fully relaxed all the atoms. 
We varied the distance between the defects and computed 
the total energies in the PM, EM and AFM states. The 
calculated exchange coupling constants are not exactly 
equal to the energy differences (SpM — -EAFM)/vacancy2^ 
since it is necessary to take into account the interaction 
of each vacancy with all periodic images. This is spe- 
cially important when the vacancy is at the edges of the 
unit cell (e.g. a vacancy located at (1,1,1), in the middle 
of the supercell, would be equidistant from 8 vacancies 
and {EpM ~ £^afm) /vacancy — 8Jiii). The resulting ex- 
change coupling constants were calculated following the 
procedure of Zhao et al.^^ and the exchange constants vs 
defects separation are shown in Fig. [J] 

We found that there are two factors that determine 
the ground state between any pair of defects, the Vsn- 
Vsn separation and the relative orientation of the anion 
adjacent to them. There is a strong competition between 
FM and AFM direct exchange interactions depending on 
the Vsn-Vsn separation when the direct coupling between 
the defects is taken into account. The total magnetic mo- 



ment of the AFM is zero for all cases except for Vsn-i (as 
marked in Fig. 2]) where we found ferrimagnetic (FIM) 
coupling. We see the defects are strongly coupled ferro- 
magnetically when the Vsn are separated by ^5.5 A. Our 
calculations also show that most of the total energy differ- 
ences between FM and AFM are very large (above RT), 
which indicate that the high Tqs in Sn02 based DMS 
come from the Sn vacancies. In Ref.'^^ the authors linked 
the observed magnetism in Sn02-based systems with the 
(P phenomenon. They extrapolated the magnetic mo- 
ments vs film thickness and assumed that if cation va- 
cancies were magnetic then the separation between them 
would be 0.5 nm, which is indeed what we observe^. This 
discloses the origin of GMM and RT ferromagnetism is 
Sn vacancies. 

Now we address the question why the distance between 
FM defects is strongly localized around 5.5 A. The Sn va- 
cancy induces a magnetization which is produced mainly 
by the O atoms. On the other hand the Sn atoms couple 
antiferromagnetically. This situation resembles a con- 
ventional superexchange mechanism, where the TM has 
some magnetic moment and the magnetism is mediated 
by the intervening O atoms. Wee see that when the an- 
gle between O atoms (0-Sn-O) is 90° the separation be- 
tween the Vsn is strongly localized around 5.5 Aand the 
most stable magnetic configuration is FM. If the angle 
between O atoms is increased to 180° then the calcula- 
tions also show FM behavior. This can not be explained 
by the typical Goodenough-Kanamori- Anderson (GKA) 
because the magnetic vacancies are coupled by more than 
one atom. Indeed, according to the GKA one would ex- 
pect the AFM behavior to be more stable for 180°, which 
is not what we find. The oscillatory behavior observed in 
Fig. m is typical of a RKKY interaction, which appears 
when the density of defects is high and the donor states 
merge with the bottom of the conduction band. By using 
the typical RKKY expressions^, 

s'm{2kFr) ~ 2kpr coa{2kFr) 



we fit our data to the RKKY model, which gives a Fermi 
wave vector kp ~ 1.19 A~^ and explains why the most 
stable FM configuration relative to the AFM configura- 
tion corresponds to a distance of 5.5 A. The oscillatory 
behavior is well defined but the amplitude of the fit is not 
well characterized due to the fact that many points are 
close to the nodes of the curve. Such oscillatory behavior 
was also shown for Co-doped Sn02 using first-principles 
calculations'^^ . A note of caution should be added how- 
ever because the high density of states at the Fermi level, 
which prompts to localized holes at the vacancies, is not 
exactly compatible with the host-like-holes limit where 
the RKKY model is supposed to be valid^. 

To further explain why the defects couple ferromagnet- 
ically, antiferromagnetically or ferrimagnetically, we take 
into account the angle between them. They may be con- 
nected directly (as shown above) or through a different 
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FIG. 5: Mechanism behind the FM and localization of the two 
defects in the Vsn system. The arrows are magnetic moments. 



path followed by different angles between them as shown 
in Fig. [3 

In Fig.[5]the Vsn marked "Fixed" remained fixed while 
1, 2, 3, and 4 are Sn vacancies at different sites. We see 
all the defects are connected through an 0-Sn-O path, i.e. 
VFixcd-0-Sn-0-V2_4, except one case where defects share 
a common oxygen atom. When we consider the coupling 
between Vpixed and V2 (VFixed-0-Sn-0-V2), for example, 
all the Sn atoms at site 1, 3, and 4 are negatively polar- 
ized and the ground state may be FM or AFM. If the 
defects share a common oxygen, i.e. Vpixod-O-Vi, then 
the ground state may be PM or FIM, depending on the 
distance between the defects. When the distance is small 
the ground state is PM. On the other hand, when we in- 
crease the distance or the angle between the defects, the 
ground state is found to be FIM, as marked (Vsn-i) in 
Fig. m The charge transfer is larger for the shared oxygen 
atom compared with other oxygens and such atom gives 
a very large magnetic moment (~ 2.00 /is) in the FM 
calculation. The AFM calculation converges to FIM and 
gives 2.00 ^s/unit cell (1.00 ^s/defect). The structural 
distortion was found to be larger for this case as com- 
pared to other cases. The angle between Vpixed-O-Vi is 
128°, which is a typical angle for which FIM has already 
been observed in spinals'^^. More recently it has been 
found theoretically that Fe-doped Sn02 also showed FIM 
(1.00 fiB /Fe) when the Fe impurities connect through O 
atoms-3>. On the other hand, if we decrease the sepa- 
ration between defects while sharing a common O atom 
then the total energy difference E-pM — Eafm becomes 
very small and PM becomes the most stable state. 

In conclusion, we used density functional theory within 
the local density approximation (LDA) and the general- 
ized gradient approximation (GGA) to understand the 
origin of giant magnetic moments (GMM) in Sn02. We 
observed that O vacancy (Vq) is non-magnetic, but Sn 
vacancy (Vsn) showed ferromagnetic behavior with a 
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large magnetic moment. Both all electron and pseudopo- 
tential approaches gave the same results. The calculated 
magnetic moment due to Vsn in Sn02 was mainly con- 
tributed by the O atoms surrounding Vsn whereas the 
Sn atoms coupled antiferromagnetically. We also calcu- 
lated the magnetic coupling between Sn vacancies and we 
showed that these defects can couple ferromagnetically, 
antiferromagnetically and ferrimagnetically. Strong fer- 
romagnetic coupling between Sn defects was found when 
they were separated by ~ 5.5 A. We also showed that the 
exchange coupling between these defects oscillates and 



such oscillatory behavior was attributed to the typical 
RKKY oscillation. The results of the extensive LDA cal- 
culations are consistent with RT ferromagnetism and the 
observed GMM, provided that an extremely high concen- 
tration of Sn vacancies of 12% (i.e. 2 vacancies per 16 Sn 
sites) would be present. 
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